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ABSTRACT 

The two-year COBE-DMR 53 and 90 GHz sky maps, in both galactic and ecliptic 
coordinates, are used to determine the normalisation of inflationary universe models 
with a flat global geometry and adiabatic density perturbations. The appropriately 
normalised cold and mixed dark matter models and cosmological constant dominated, 
cold dark matter models, computed for a range of values of Qf, and h, are then com- 
pared to various measures of structure in the universe. Critical density CDM models 
appear to be irreconcilable with observations on both large and small scales simultane- 
ously, whereas MDM models provide a somewhat better fit to the data. Although the 
COBE-DMR data alone prefer a nearly critical value for the total density, low-density 
cosmological constant models with fi > 0.15 can not be rejected at a confidence level 
exceeding 95%. Such models may also provide a significantly better fit to the matter 
distribution data than critical density CDM. 
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The COBE-DMR discovery of cosmic microwave hack- 
le ground (CMB) anisotropy (Smoot et al. 1992, Bennett et 
1992, Wright et al. 1992) has affected cosmology in 
Q both epistemological and practical ways, but its predomi- 
nant quantitative influence has been to provide the means 
t/3 for the accurate normalisation of theories of large scale stree- 
ts ture formation. 

The development of inflationary ideas during the 1980s 
induced over a decade-long adherence to the cosmological 
paradigm which posits that the universe is spatially flat. 
Such a picture requires that the present energy density of 
the universe is dominated by non-baryonic dark matter or 
alternatively by a non-zero vacuum energy contribution (a 
cosmological constant term, CIa). 

The minimal version of the model, which invokes cold 
dark matter (CDM) as the major constituent of the uni- 
verse, is presently in direct confrontation with astronom- 
ical observations. Extensions of the model, which in ad- 
dition to CDM postulate either an admixture of hot dark 
matter (HDM) or a non-zero cosmological constant (CDM- 
A), enjoy considerable popularity in contemporary cosmo- 
logical research. Whilst vigorous discussion ensues in the 
literature as to the plausibility of the mixed dark matter 
(MDM) (see e.g. Schaefer, Schafi & Stecker 1989, Klypin et 
al. 1993, Pogosyan & Starobinsky 1993, Ma & Bertschinger 
1994, Primack et al. 1995) and CDM-A (e.g. Peebles 1984; 
Efstathiou, Maddox, & Sutherland 1990: Carrol, Press, & 
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Turner 1992; Kofman, Gnedin, & Bahcall 1993; Krauss & 
Turner 1995; Ostriker & Steinhardt 1995) models as viable 
cosmologies, we consider both CDM and MDM versions of 
the inflationary scenario (with Qa = 0) and CDM models 
with a non-zero cosmological constant. 

In this paper, we use the linear angular power spectrum 
estimation technique (Gorski 1994) to normalise dark mat- 
ter models to the two year COBE-DMR anisotropy mea- 
surements (Bennett et al. 1994). Bunn, Scott & White 
(1995), Bunn & Sugiyama (1994) and Sugiyama (1994) have 
used different power spectrum normalisation schemes to ob- 
tain comparable results for those models in common with 
this paper. We then discuss some predictions for large-scale 
structure measures resulting from this normalisation. 



2 CMB ANISOTROPY NORMALISATION 
PROCEDURE 

2.1 Data Selection and Power Spectrum Inference 
Method 

The COBE-DMR project generates sky maps from the 
time-ordered data in both galactic and ecliptic reference 
frames. Previously, published work by the DMR team has 
utilised only the galactic data sets, whilst the general com- 
munity analysed the publically released sky maps in ecliptic 
co-ordinates. (Note that both sets of data are now avail- 
able by anonymous ftp from the NSSDC). In this work, the 
two year COBE-DMR 53 and 90 GHz galactic sky maps 
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Figure 1. CMB anisotropy power spectra normalised to COBE-DMR, where A| = 1(21 + l)a|/47r. a) CDM/MDM: the four sets of 
curves correspond to the power spectrum amplitude determination as follows (from top to bottom): 1) ecliptic coordinates, quadrupole ex- 
cluded, 2) ecliptic coordinates, quadrupole included, 3) galactic coordinates, quadrupole excluded, and 4) galactic coordinates, quadrupole 
included. Within each set the three curves denote: h — 0.5 and f2j = 0.05 - heavy type, h — 0.3 and f2(, = 0.14 - medium-heavy type, 
h = 0.8 and Qi, = 0.02 - light type, which are consistent with the constraints from BBN. The crossing point within each set is at I tss 11. 
Several pure Sachs- Wolfe power law spectra are shown (shifted down to 0.9 at I — 2). Although a power law approximation, P(k) oc k n , 
used to generate the multipole coefficients a 2 , solely through the Sachs & Wolfe (1967) effect will be poor for such spectra, a value of 
n ~ 1.1 would be most appropriate over the range I ;J 15. This is a little steeper than the underlying, inflationary n — 1 spectrum (Bond 
1993). The proposed overall two year COBE-DMR normalisation (translated from Q rma -ps ~ 20 iiK) is represented by a filled circle; 
the error bar represents a typical statistical uncertainty of an individual likelihood fit. Clearly, this encompasses the fit uncertainties due 
to choice of coordinate system, quadrupole inclusion/exclusion, and/or cosmological parameter values. 

b) CDM-A : the thin lines correspond to the h — 0.5, f2j = 0.05 models, thick lines to models with h — 0.8, f2(, = 0.02. The curves are 
drawn for four values of total density (in descending order at I = 2): Qq = 0.1, 0.2, 0.3, 1. The approximate pivot point (for Qgh 2 Jj 0.05 
models) is also shown. 



are used identically as in Gorski et al. (1994). In addition, 
we also utilise the ecliptic coordinate frame data sets as a 
check on the extent to which the coordinate dependent noise 
binning can affect the results of the analysis. 

We implement the power spectrum estimation method 
as described in considerable detail in Gorski (1994). Co- 
ordinate system specific orthogonal basis functions for the 
Fourier decomposition of the sky maps are constructed so as 
to include exactly both pixelisation effects and a |6| < 20° 
galactic plane excision (leaving 4016 and 4038 pixels in the 
galactic and ecliptic sky maps, respectively). A likelihood 
analysis is performed as in Gorski et al. (1994). 

2.2 Theoretical Spectra of Anisotropy 

We consider spatially flat models with fio + Oa = 1 and pri- 
mordial perturbations as prescribed (to within an arbitrary 
amplitude) by the inflationary scenario. Hence, random- 
phase, Gaussian, scalar curvature perturbations are assumed 
(no gravity waves), corresponding to adiabatic density per- 
turbations with the Harrison-Zel'dovich power spectrum, 
P(k) oc k. 

For fio = 1, the bulk mass density is provided by either 
cold or mixed dark matter; in the case of the MDM models, 
the hot dark matter is introduced in the form of either one 
or two (equal mass) families of massive neutrinos, with a 
contributed fraction of the critical density Q„ = 0.15, 0.2, 



0.25 and 0.3 for one massive flavour, and Q. v = 0.2, 0.3, 
otherwise. For the CDM model, the values of the Hubble 
constant, H = 100 h km s _1 Mpc -1 , and baryon abun- 
dance are sampled at h = 0.3, 0.4, 0.5, 0.6, 0.7, and 0.8, 
and fi 6 = 0.01, 0.03, 0.05, 0.07, 0.1 (for all values of h), re- 
spectively; in order to trace the Big-Bang nucleosynthesis 
(BBN) relation, Q b = 0.013 /i" 2 (Reeves 1994), we also use 
Q b = 0.02 for h = 0.8, and Q b = 0.14 for h = 0.3. For the 
MDM models, we choose h = 0.5 and 6 = 0.05. The CDM- 
A models considered are specified by either h = 0.5 or 0.8, 
and Q b is again selected to match BBN constraints, i.e. the 
baryon density is Q& = 0.05 for h=0.5, and 06 = 0.02 for 
h = 0.8. For comparison we also consider CDM-A models 
with fif, = 0.03 for both values of h. 

The CMB anisotropy multipole coefficients and the 
matter perturbation transfer functions for both models were 
evaluated using the Boltzmann equation code of Stompor 
(1994) by solving the propagation equations up until the 
redshift z = 0. Thus all effects caused by the cosmological 
constant were taken into account exactly. 

Fig. la shows a selection of the CDM radiation power 
spectra normalised to the COBE-DMR two year CMB 
anisotropy. Over the low-^ range of CMB multipoles probed 
by COBE-DMR the theoretical spectra are indistinguish- 
able for fio = 1 CDM and MDM models with equivalent h 
and fif,. In this case, the power spectrum amplitude derived 
from the data applies equally to both CDM and MDM mod- 
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Figure 2. a) Contour plots of the likelihood density in the Q r ms-ps — &o plane for the h = 0.8 CDM-A models determined from 
analysis of the two year COBE-DMR data in ecliptic coordinates. 68%, 95%, and 99.8% integrated likelihood levels, derived under the 
assumption of a uniform (Q r ms-PS, f2o)~prior, and the Qq =const-ridge lines are shown. Thick lines correspond to the results of the 
analysis with the quadrupole included, thin lines excluding the quadrupole. Diamonds show the loci of the likelihood maxima: filled 
symbols show the situation when the quadrupole is included, open symbols when it is excluded. The filled and open circles show the 
likelihood maxima loci derived from analysis of the galactic data to illustrate the small shift in the inferred model amplitude. 

b) As in a for h = 0.5. 

c) Marginal likelihood densities for Qq, derived by integration of the two dimensional density distributions (from panels a and b) over 
Qrms-ps, are shown by thin lines. The functions corresponding to the analysis with the quadrupole included fall more steeply with 
decreasing f2o than when the quadrupole is excluded. The lines for the h = 0.8 model extend down to f2o = 0.02. Results derived with 
both the ecliptic and galactic frame data are shown — the ecliptic frame results fall below the galactic ones at small f2o- The heavy lines 
show the cumulative likelihood functions derived from the above likelihood densities. The quadrupole included case renders more steeply 
rising curves with h = 0.5, 0.8 models, and galactic/ecliptic results all merge together. The less steeply rising curves, split at very small 
values of f2o, correspond to the no quadrupole case. The 5% cumulative likelihood level is shown, as is the overall most restrictive, h- 
and coordinate frame-independent, Qq > 0.15 limit on the CDM-A models implied by the two year COBE-DMR data alone. 



els. Fig. lb corresponds to several CDM-A radiation power 
spectra. At \ow-£ 30), the radiation power spectra are 
determined by the usual Sachs & Wolfe (1967, SW) effect 
and by two integrated SW effects (ISW), a low-redshift one 
- induced by the non-zero cosmological constant, (Kofman 
& Starobinsky 1985, Gorski, Silk, & Vittorio 1992), and a 
high-redshift one - due to the non-negligible radiation en- 
ergy density contribution to the total density after decou- 
pling. For those fc-modes which dominate the ai coefficients 
with £ & 10, the high-redshift ISW effect simply cancels 
that part of the SW effect due to radiation density pertur- 
bations present on the last scattering surface. The sum of 
both effects therefore scales with cosmological parameters as 
would a matter perturbation generated potential. In partic- 
ular, its dependence on h is almost the same as that of the 
low-redshift ISW contributions, while the dependence on fio 
differs because of the fio-dependence of the density pertur- 
bation growth rate at low redshifts. Hence, the shape of that 
part of the radiation transfer function depends very weakly 
on the Hubble constant. However, this is not the case over 
the ^-range 10 £ 30, where the canceling of the high-z 
ISW and "radiation" part of the SW effect is less complete 
and the dependence of the power spectrum on h is clearly 
noticeable. The strength of both ISW effects falls with in- 
creasing fio- Consequently, with respect to the fiducial pure 
Harrison-Zel'dovich form ~ !//(/+!), the power spectra fiat- 



ten as fio increases, and for models with near critical matter 
density the spectra monotonically grow with £. This steep- 
ening of the low-^ end of the power spectrum over the pure 
SW case is caused by the high-redshift ISW effect, which, 
though small, does not vanish even in critical density cases 
(Bond 1993; Bunn et al. 1995), and a small contribution 
from the first "Doppler" peak. The steepest spectra arise in 
models near critical matter density (Oo — 0.8). 



We parameterise all such spectra by the exact value 
of the quadrupole, Qrms-ps (the value of the exactly 
computed quadrupole, 02, multiplied by (5/47r) 1 ^ 2 , which 
depends non-trivially on cosmological parameters), in a 
straightforward generalization of the Q rm s-ps introduced in 
Smoot et al. (1992) for pure power law model spectra. When 
Ho ^ 0.5, the ISW contributions boost the quadrupole over 
its pure Sachs & Wolfe counterpart. For fio ^ 0.5 — models 
with a negligible cosmological constant induced ISW effect 
— the resulting exact quadrupole is lower than its pure SW 
value (see Table 2). In the case of critical density models, 
these differences, mainly due to the high-redshift, integrated 
SW effect, are rather small, being equal to ~ 10%, 5%, or 
2% for h ^ 0.4, h ~ 0.5, and h ~ 0.8, respectively, of the 
pure SW value. 
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2.3 Results of Qrms-ps Fitting 

A typical likelihood fit of flat dark matter models to the 
two year COBE-DMR data yields a ~ 13a significant de- 
termination of Qrms-ps, with a value of ~ 20 fiK for the 
critical matter density case, and in the range ~ 20 — 26 /iK 
for Qa > 0. Systematic shifts in the central value of the 
fit are observed due to: 1) differences in the noise pixeli- 
sation in the galactic and ecliptic coordinate frames, which 
result in a ~ 0.8 fiK difference between the inferred normal- 
isation amplitudes, with higher values obtained from the 
ecliptic maps (an effect alluded to, but not quantified, by 
Bunn et al. 1995); 2) exclusion of the quadrupole, which 
produces a ~ +0.4 fiK variation in the fitted amplitude; 3) 
the uncertainty in the values of h and 06 (reflected by differ- 
ences in spectral shape over the Grange accessible to DMR) 
causes an additional small spread, ±0.25 /xK, of the fitted 
amplitudes. 

Likelihood contours in the Qrms-ps — Oo plane (under 
the assumption of a uniform prior) are shown in Figs. 2a, 
b. The preferred value of the matter density is clearly close 
to critical, but even low-density models with fio > 0.15 lie 
well within the 95% confidence limits and appear acceptable. 
Due to the weak dependence of the COBE-DMR normalised 
Qrms-ps on h, this 95% C.L. lower limit is essentially Hub- 
ble constant independent. 

A convenient summation of the proposed overall nor- 
malisations for the flat, dark matter models is then: 

Qrms-Ps{liK) ~ [20 + g(Q ,h)\ x (1± 0.077) 
±0.4 ±0.2 ±0.25 

where g{Q ,h) = 9.1/i 06 exp(-O§/i - 3 /0.057) for O < 1, 
and <?(fio, h) = for fio = 1. The error ranges represent the 
statistical error and uncertainties associated with effects 1 
through 3 above, respectively. It will be noted that the sta- 
tistical error on the inferred normalisation is considerably 
larger than the other uncertainties. Exact values for a rep- 
resentative set of CDM-A models can be found in Table 2, 
which can be easily interpolated to other models with high 
accuracy. However, our fit is good to within 5%. 

Gorski et al. (1994) showed that for the power law 
models specified by Qrms-ps and the spectral index n a 
convenient, n-independent normalisation was in terms of the 
multipole amplitude ag ~ 8 fiK. This border is related to 
the point at which the theoretical signal to noise ratio is 
of order unity. For the O = 1 dark matter models, we 
find that the appropriate pivot point is at an ~ 7.15 ± 
0.55/iK. However, those models with a non-zero cosmological 
constant show no pivot point at a common multipole order. 
Nevertheless, for those models with Qoh 2 <; 0.05, a good 
approximate pivot can be found at I ~ 8 (Fig. lb), yielding 
the model-independent normalisation as — 9.5 ± 0.73/xK. 

3 MEASURES OF LARGE-SCALE 
STRUCTURE 

Having determined the normalisations for our grid of 
models, we now proceed to discuss the values of several 
large-scale structure statistics. These have been computed 
from the matter perturbation spectra according to the usual 
prescriptions given as footnotes to Table 1 (which contains 
representative values for the fio = 1, h = 0.5, fi& = 0.05 



Table 1. Inferred cosmological statistics for models with f2o = 
V-cdm + Sl v +Q b = l,n b = 0.05, h = 0.5, and a COBE-DMR 
normalisation of Q r ms-PS = 20.04 jitK. N„ is the number of 
massive neutrino species, and m„ the neutrino mass in eV. The 
errors, including both statistical (1<t) and systematic deviations, 
are of the order of 11%. 
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models). Table 2 summarises the pertinent values for a set 
of cosmological constant models. 

3.1 Mass Fluctuations: as and J3 

COBE-DMR normalised values of (as)mass are shown 
in Fig. 3a. A related, observable quantity is the fluctuation 
in the number of galaxies within a sphere of fixed radius. 
Recent estimates are close to the standard value from Davis 
& Peebles (1983), (<Jg) ga i — 1 — a representative selection of 
the galaxy surveys is: (as) ga i — 0.8 (Fisher et al. 1993), 0.87 
± 0.07 (Feldman, Kaiser & Peacock, 1994), and 0.83to;o7 
(Baugh & Efstathiou, 1993). 

In order to compare the rms mass fluctuations with the 
rms galaxy fluctuations, we adopt the simple picture of bi- 
ased galaxy formation, wherein (as) ga i = b (<T8)mass, and 6 is 
the linear bias factor. Some constraints on the value of 6 can 
be imposed by recent galaxy surveys, but are more uncer- 
tain than the (as) g ai determinations and usually inseparable 
from the total density parameter. 

In the fio = 1 case, the estimated values of 6 range from 
0.9 to 2, with strong bimodal behaviour about 6 ~ 1 and 2 
(see Table 1 in Dekel 1994, and Table 1 of Cole, Fisher & 
Weinberg 1994, with fio = 1). Estimates of (as)mass in- 
ferred from galaxy cluster catalogues favour a higher value 
for 6. Henry & Arnaud (1991) used the abundance of clus- 
ters as a function of X-ray temperature to derive (as)mass = 
0.59 ± 0.02 for a scale-free, flat universe. White, Efstathiou 
& Frenk (1993) have used the masses and abundances of rich 
clusters of galaxies to determine (as)mass — 0.52 - 0.62 for 
a critical density universe (relatively independent of ft). A 
low value of (as)mass ^ 0.5 is also required in critical den- 
sity models by the observed low pair-wise velocities (Davis 
et al. 1985, but see Zurek et al. 1994). On the other hand, 
Seljak & Bertschinger (1994) determine a 95% confidence 
range for (a 8 ) m£ISS of 0.7 - 2.3 from the POTENT density 
field reconstruction. 
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Table 2. Inferred cosmological statistics for chosen spatially flat CDM-A models with 
Qb = 0.013h -2 , and a COBE-DMR normalisation, expressed in terms Qrms-PS- The 
errors, including both statistical (la) and systematic deviations, are of the order of 11%. 
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a 2(sw) is the pure Sachs- Wolfe contribution to the quadrupole. 



Although the value of the rms mass fluctuations is not 
yet well known, it seems unlikely that the biasing parameter 
is less than unity. The generally high values of (as)mass > 1 
predicted from the COBE-DMR normalised O = 1, CDM 
models are difficult to reconcile with observations. Physical 
motivations for requiring a bias parameter somewhat larger 
than unity were pointed out early on (e.g. Silk 1977), and 
are supported by recent hydrodynamical simulations (Katz, 
Hernquist & Weinberg 1992, Cen & Ostriker 1992) which 
suggest a value 6 ~ 1.5. The values of 6 > 1 would require 
either pure CDM in an fio = 1 universe with a low-/i or 
high-fij, in disagreement with observations (0.4 & h ^ 0.8) 
or the BBN constraints, respectively, or alternatively an ap- 
preciable admixture of a hot dark component. 

A non-vanishing cosmological constant can also pro- 
vide a means to circumvent these problems. The predicted 
{o's)mass as a function of Qoh 2 is shown in Fig. 36. A low 
value of fio suppresses peculiar motions so high bias is not 
required to fit the small-scale galaxy pair-wise velocities. 
Thus unbiased models with (as)mass — (o's)gai are likely to 
be of interest. 

Note that a high biasing parameter (6 J> 1.5) to- 
gether with the (a s )mass ^ 0.57Q 66 (White et al. 1993, 
Strauss & Willick 1995), estimation based on galaxy clus- 
ter catologs would require an almost critical density model 
with a low (h ^ 0.3) Hubble constant in order to be com- 
patible with the COBE-DMR measurement. The situa- 
tion is even more dramatic for such models if the Seljak 
& Bertschinger (1994) estimation of rms mass fluctuations, 
(o's)mass — 1 .3 ^J'g Qq 6 , is considered. In this case, none of 
the highly biased (6 ^ 1.5) subcritical total density models 
can fulfill this constraint, and even the relaxed requirement 
that 6 > 1 imposes a 95% C.L. bottom limit on the total den- 
sity Oo ^ 0.55 which, in conjunction with the COBE-DMR 



data, again imposes an upper limit on h (^ 0.5). Until the 
apparent discrepancy between the two (as)mass estimates is 
resolved, a conservative approach is adopted with 6 > 1 as 
suggested above, so that models with Q h 2 0.13 are also 
relevant. 

Comparison of the J3-integral, computed on the larger 
scale of R = 20/i _1 Mpc (which, in principle, is less con- 
taminated by non-linear galaxy power spectrum evolution 
than as), to the observational value from Davis & Peebles 
(1983), J 3 (20/i _1 Mpc) ~ 700 h~ 3 Mpc 3 (probably accurate 
to within ~ 30%), also seems to favour the low-fio CDM-A 
-models (Fig. 4) in agreement with the above constraint. 

3.2 Large-Scale Flows 

The local streaming motions of galaxies provide an interest- 
ing constraint for cosmological models. In particular, galaxy 
peculiar velocities directly measure mass fluctuations, inde- 
pendently of a linear bias parameter. Dekel (1994) gives 
estimates of the average peculiar velocities within spheres 
of radius 1000 to 6000 km s -1 , which seem to be in agree- 
ment with the HI data from Giovanelli & Haynes (Dekel, 
private communication). The COBE-DMR normalised rms 
bulk flows are shown in Fig. 5. The critical density model 
predictions are in good agreement with the POTENT data, 
and it is apparent that all models with Qoh ^ 0.06 can 
account for the POTENT velocities. However, an impor- 
tant recent development in the field was the determination 
by Lauer & Postman (1994) of significant bulk flow in a 
deep volume limited sample of ~ 100 galaxy clusters. If 
confirmed, even the critical density models will be in serious 
trouble (since all of the presently considered models of struc- 
ture formation predict too rapid a decrease in the bulk flow 
amplitude with scale. See also Feldman & Watkins 1994, 
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Figure 3. a) (ag) ma ss values predicted from the COBE-T)MR normalised, flat, dark matter dominated models. The thin solid lines 
correspond to the CDM model with h decreasing from 0.8 to 0.3 in steps of 0.1 from top to bottom. The thick solid line shows the 
models which obey the Big Bang nucleosynthesis constraint, f2(, = 0.013 /i -2 . The individual points represent several mixed dark matter 
models: triangle - Q. v — 0.15, square - Q. v — 0.20, hexagon - f2„ — 0.25, diamond - Q. v — 0.30, where the filled symbols have N v = 1, 
open symbols N v = 2. All of the points correspond to h = 0.5 and f2j = 0.05 but are spread out on the plot for clarity. 
b) (&8)mass values predicted from the COBE-DMR normalised, CDM-A models, plotted as a function of Qgh 2 . The h = 0.5, f2(, = 0.05 
models are shown by solid lines terminating at Qoh 2 = 0.25, while the h = 0.8, f2j = 0.02 lines extend to the edge of the plot. For 
each model, the upper of the pair of lines corresponds to the ecliptic frame/no quadrupole analysis, whilst the lower represents the 
galactic frame/quadrupole included case (the remaining two combinations would lie in between). Dotted lines show the f2j = 0.03 results 
averaged over the four possible combinations. 




Figure 4. a) The values for the J3 integral over the density perturbation correlation function within 20 h -1 Mpc (units [h -1 Mpc] 3 ). 
Same coding as 3a. 

b) The values for the J3 integral over the density perturbation correlation function within 20 fe -1 Mpc (units [h -1 Mpc] 3 ). Same coding 
as 36. 



Strauss et al. 1995.). 

3.3 Galaxy and Matter Distribution Power 
Spectra 

In principle, the most direct and informative compari- 



son of theory and observations can be conducted using the 
predicted and empirical power spectra. Several recently es- 
timated galaxy power spectra are shown in Fig. 6, together 
with various theoretical linear power spectra (assuming 6 
=1). However, such a comparison is still subject to difficul- 
ties, including the linear nature of our theoretical computa- 
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— top, h = 0.5, f2(, = 0.05 — middle, h = 0.3, f2(, = 0.14 — bottom. Thin lines separating from the h = 0.5 curve correspond to the 
MDM models with N„ = 1 (higher line), and iV„ = 2, practically independent of f2„. 

POTENT data (courtesy A. Dekel) are shown by circles, and the square shows the Lauer & Postman datum. 

b) rms amplitudes of the large scale flows ([kms -1 ]) as a function of the top-hat sphere radius encompassing the sample (an additional 
r s = 1200 km s _1 smoothing is included). Heavy lines — h = 0.8, Q b = 0.02, and Q = 1, 0.4, 0.3, 0.2, 0.1, 0.05, and 0.02 in descending 
order. Thin lines — h = 0.5, Q b = 0.05, and do = 1, 0.4, 0.3, 0.2, 0.1, and 0.05 in descending order. 

c) rms amplitudes of the bulk flow ([kms -1 ]) within spheres of 40, 60, and 100 h -1 Mpc plotted in descending order, respectively, as a 
function of Ugh. Heavy lines — h = 0.8, f2o = 0.02, thin lines — h = 0.5, f2o = 0.05. 



tions (since on small scales, the effects of non-linear evolu- 
tion should be apparent), possible distortions introduced by 
redshift-real space mapping, and the lack of agreement be- 
tween different observational estimates of the galaxy power 
spectrum (which may be explicable in terms of intrinsic vari- 
ations in the survey samples and/or by differences in the 
estimation technique utilised), fn spite of such difficulties, 
we still consider such a comparison to be informative in a 
qualitative way. 

It is apparent from Fig. 6a that none of the empirical 
real space power spectra (Baugh & Efstathiou 1993, Peack- 
ock & Dodds 1994) exhibit as much small-scale power as 
the linear COBE-VMR normalised fi = 1 CDM spec- 
trum. (Note, however, that the redshift space spectra 
from da Costa et al. (1994) indicate more power on small 
scales than the real space APM spectrum, a discrepancy 
which can only become more pronounced when one recalls 
that the real-to-redshift space mapping suppresses power on 
scales k <; 0.1/iMpc -1 - Gramann, Cen & Bahcall 1993 - 
and would require high, scale dependent bias.) This could 
be considered as supportive of the MDM models. How- 
ever, the high amplitude of the matter power spectra im- 
plied by the COBE-DMR measurements leaves very little 
room for the biasing parameter to exceed unity, in conflict 
with the 6-estimates from cluster properties. Baugh & Efs- 
tathiou (1993) have suggested that, with b ~ 2, the matter 
power spectrum required by the observed galaxy distribu- 
tion should look rather like that for the fio = 1, T = 0.2 
model (here, we follow the parametrization of the matter 
power spectra introduced by Efstathiou, Bond, & White 
1992) depicted in Fig. 6a as a thin solid line. Alas, none of 
the COBS-DMR-normalised CDM or MDM model power 



spectra resembles such an ad hoc spectrum. 

Considering the non-zero A models shown in Fig. 6b, 
we see that on small scales (A: > 0.1 /iMpc -1 ), the estimated 
power spectra from Baugh & Efstathiou (1994), Peacock & 
Dodds (1994) and the IRAS 1.2 Jy redshift survey are in rea- 
sonable agreement for those models with Qoh 2 0.08. This 
limit is somewhat more restrictive than that obtained previ- 
ously, but, combined with the COBE-DMR plus POTENT 
inferred lower limits on fio, still leaves room for acceptable 
models. (If the more stringent gravitational lensing con- 
straints on Oa < 0.7 - Maoz & Rix 1993 - are taken into 
account, this range shrinks somewhat). However, since ap- 
parent discrepancies exist between different empirical power 
spectra on these scales which can not be solely ascribed to 
redshift space induced distortions, relaxation of the above 
constraint would seem to be in order until the observational 
situation is clarified. Conservatively, we retain the previous 
limit. 

However, it then appears that the preferred CDM-A 
models can not simultaneously fit the large-scale behaviour 
indicated by the APM data. This is more properly de- 
scribed by an unbiased CDM h = 0.5 model, which, un- 
fortunately, fails noticeably on smaller scales. One should 
note, however, that on these larger scales there is again pro- 
nounced disagreement between different data sets (with the 
SSRS2+CfA2 data of da Costa et al. , 1994, clearly over- 
shooting all other estimates by a factor ~ 4, which can- 
not be completely explained by the effect of projecting the 
power spectrum to redshift space. Hence, it would require 
b <> 2fi[!- 66 /5 < 1 - Kaiser 1987. For values of 6 in the 
range ~ 1 — 2, the linear correction appropriate for scales 
k <, 0.03 AMpc" 1 - Gramann et al. 1993 - is only 10%-30% 
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Figure 6a. COBE-DMR normalised inhomogeneity power spectra and, shown as the lowermost curves, miscellaneous spectral windows 
required for the computation of the statistics considered in this paper (as annotated). Overplotted are several observational estimates 
for the galaxy distribution power spectrum: squares - Baugh & Efstathiou (1993), filled circles - Peacock & Dodds (1994), pentagons - 
Fisher et al. (1993), open and filled triangles - two data subsets from da Costa et al. (1994). The theoretical mass distribution power 
spectra for the Qq = 1, h = 0.5, f2(, = 0.05 models are shown in the middle section of the plot: the top heavy line shows the CDM 
power spectrum; the lower heavy lines correspond to two N v = 1 MDM models — Q v = 0.2 (upper) and f2„ = 0.3 (lower); the medium 
heavy lines correspond to the equivalent N„ — 2 MDM models. Note that the MDM transfer function for N v — 2 drops faster near 
k ~ 0.04 /iMpc -1 N„ = 1. The thin line shows the T = 0.2 transfer function for CDM. The vertical error bar on the P(k) ~ k part of the 
spectrum above the an window illustrates the ~ 1 a allowed variation in amplitude for this fixed slope. Conversely, the superposed 'bow' 
attempts to indicate the allowed 1 cr-variation (~ 0.3) in tilt for spectra with the fixed an = 7.15 iiK amplitude. Unfortunately, there is 
no obvious way to represent the two dimensional uncertainty in the spectrum determination from the two year COBE-DMR data for 
flat dark matter models. One is not at liberty to arbitrarily vary both amplitude and spectral slope simultaneously. This is equivalent 
to the observed degeneracy seen in the two-dimensional (Q r ms — PS> n ) nls t0 pure power law models (Gorski et al. 1994 and references 
therein). The upper section of the plot reproduces the theoretical spectra for CDM, MDM (N„ = 1, Q„ = 0.2), and T = 0.2-CDM - all 
shifted upward for clarity. 



for fi = 0.1 and 20%-50% for fi = 0.4 so the pronounced 
power enhancement survives.). It undoubtedly remains too 
early for any definitive conclusions to be reached. If the 
SSRS2+CfA2 data are confirmed, the power enhancement 
visible on scales k ~ 0.01 Mpc -1 may favour the low-fio 
CDM-A models. Improved, very deep galaxy catalogs may 
provide us with a definitive test of CDM-A models. 



4 DISCUSSION 

The improved quality of the two year COBE-DMR data 
combined with reliable power spectrum estimation tech- 
niques allows the accurate normalisation of cosmological the- 
ories. Previously, the variance of COBE-DMR temperature 



fluctuations on a 10° angular scale was utilised for the nor- 
malisation of the power spectrum. Subsequent work (Wright 
et al. 1994, Banday et al. 1994) has demonstrated that 
this technique can be unreliable without considerable atten- 
tion. More appropriate methods take full advantage of the 
measured CMB anisotropy power distribution on all angular 
scales accessible to the COBE-DMR instrument, as imple- 
mented in this paper. 

The cold dark matter theory with a standard choice of 
cosmological parameters requires a very high normalisation 
in order to fit the CMB anisotropy distribution. Analy- 
sis of the first year of COBE-DMR data had already sug- 
gested that (tJ8)mass ~ 1 (Wright et al. 1992, Efstathiou 
et al. 1992), and this value increases to ~ 1.4 with two years 
of data and an improved analysis technique. Although this 
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Figure 6b Same as figure 6a, but for non-zero cosmological constant models. Solid lines, in descending order, show the CDM-A power 
spectra for n = 0.1, 0.2, 0.3 — heavy type, and Q = 0.1, 0.2, 0.3, and 0.4 — light type. The h = 0.5, Q b = 0.05, and h = 0.8, Q b = 0.02 
models are shown in thin and heavy lines, respectively. Dotted lines correspond to f2o = 1, CDM models. Since the dependence of the 
fitted Qrms-ps amplitude on model parameters is rather weak, and because the SW and ISW effects scale with the total density, there 
is a strong resulting variation with f2o of the amplitudes of the matter power spectra, as shown in Fig. 4. Note that the quadrupole 
window function depends (due to the SW and ISW effects) on model parameters, Q 2 , ma _ PS = (27r 2 ) -1 J dk P(k)WQ(k). For the pure 
Sachs- Wolfe effect, Wq(&) = Q^H^j^ikRLSs)/^ 2 , where Rlss is the radius of the last scattering surface (at a redshift zlss — 1100). 
The windows for Qq = 0.1, 0.3, and 1 are shown from left to right (note the trend in the window maxima and changes in the shape of 
the functions). 



normalisation allows the theory to predict large-scale veloc- 
ities of comfortably high amplitude, it also results in a sig- 
nificant overproduction of density perturbations on scales 
<, 20/i _1 Mpc. CDM has often been criticised for its poor 
match to both galaxy and cluster distributions. Mixed dark 
matter models manage to circumvent, to a certain degree, 
these same problems by construction — massive neutri- 
nos partially damp the perturbations at those length-scales 
where CDM looks problematic. Among the MDM mod- 
els those with two species of massive neutrino seem to meet 
the observational constraints more comfortably (see also Pri- 
mack et al. 1995). The larger free-streaming radius allows 
for the suppression of the perturbation amplitude on larger 
scales than in models with one massive flavour. This is re- 
flected in the decrease of the predicted (<Tg)mase values. Nev- 
ertheless, the proponents of MDM will have to address the 
viability of the model viz. the simultaneous requirements 
that there should be no bias between the galaxy and mass 
distribution (as suggested by this analysis) and the galaxy 
pair-wise velocities should be small. 

The low-Oo CDM-A (with O > 0.15) models can not 



be challenged solely on the basis of the COBE-DMR data, 
although the most likely fio value is strikingly close to unity 
(see also Bunn & Sugiyama 1994). Conversely, observa- 
tions of the matter distribution do require low values of the 
total density resulting in the constraint Qoh 2 ^ 0.13. On 
larger scales, these models predict bulk flows in reasonable 
agreement with POTENT for Qoh <; 0.06, but dramatically 
smaller than determined by the Lauer & Postman analysis 
for any (< 1) value of Oo- Therefore, if the latter observa- 
tion is confirmed all CDM-A models (together with critical 
density CDM or MDM models) will be found wanting. The 
predicted mass distribution on scales of a few hundred Mpc 
overshoots most of the experimental estimates (including 
those derived from the IRAS 1.2Jy & APM surveys), but is 
in rough agreement with that recently obtained by da Costa 
et al. (1994). 

Whilst it would be safer to await the final 4-year COBE 
results before offering definitive statements as to the viabil- 
ity of theoretical models, one should note that the CDM 
normalisation derived from the two year COBE-DMR data 
does appear to be irreconcilably high, while the MDM model 
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has little room left for adjustment. Although the position of 
the CDM-A models is more comfortable than that of critical 
density CDM models, it does not remain free from poten- 
tially fatal flaws. Unfortunately, CMB anisotropics on a one 
degree scale do not offer any serious prospects for distin- 
guishing between cosmological constant models and other 
viable scenarios (Bond et al. 1994). 

Krauss & Turner (1995) have recently proposed that 
cosmological constant CDM models with Oa in the range 0.6 
- 0.7 can account for both large-scale structure and CMB 
observations. Similarly, Ostriker & Steinhardt (1995) have 
suggested that a spatially flat model with Oa = 0.65 ±0.1 
and a small tilt, 0.8 < n < 1.2, is consistent with a num- 
ber of cosmological constraints (including direct measures 
of H , stellar ages, large-scale structure statistics and CMB 
anisotropies). Accepting their challenge now to identify 
a problem for such models, we refer the reader again to 
fig. 6b. It appears that agreement may be achieved on both 
COBE and as scales. However, because the 6 > 1 condi- 
tion requires a relatively low value of the Hubble constant, 
h <, (0.13/O ) 1/2 ~ 0.65 for O > 0.3 in the Harrison- 
Zel'dovich small negative tilt (n < 1) may be nec- 

essary if h is to be as high as ~ 0.7 — 0.8. Such a tilt 
would enhance the peak in the matter power spectrum at 
k ~ 0.01 /iMpc -1 - where the current observational status 
remains uncertain - relative to the power on as scales near 
k ~ 0.1/iMpc _1 . Therefore, definite conclusions will have 
to await more reliable large-scale, deep galaxy surveys ca- 
pable of unraveling the shape of the galaxy power spectrum 
on scales up to comoving length-scales ~ 600 /i _1 Mpc. 
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